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Abstract. Perfect imaging for electromagnetic waves using the Maxwell Fish Eye 
(MFE) requires a new concept: the perfect drain. From the mathematical point of view, 
a perfect point drain is just like an ideal point source, except that it drains power from 
the electromagnetic field instead of generating it. We show here that the perfect drain 
for the MFE can be seen as a dissipative region the diameter of which tends to zero. The 
complex permittivity s of this region cannot take arbitrary values, however, since it 
depends on the size of the drain as well as on the frequency. This interpretation of the 
perfect drain connects well with central concepts of electromagnetic theory. This opens 
up both the modeling in computer simulations and the experimental verification of the 
perfect drain. 



1. Introduction 

The possibility of focusing light below the diffraction limit (super-resolution) has been 
demonstrated in the last decade using left-handle materials [1][2][3][4] (that is, materials with 
negative dielectric and magnetic constants [5]). Recently, a new possibility has been introduced 
using the Maxwell Fish Eye (MFE) lens, for a material with a positive, isotropic and 
inhomogeneous refractive index. It is well known that, in the Geometrical Optics (GO) 
framework, the MFE perfectly focuses rays emitted by an arbitrary point of space into another 
(its image point). Leonhardt [6] has demonstrated in that the MFE lens in two dimensions not 
only perfectly focuses radiation in the GO approximation, but also does so for actual fields of 
any frequency, a result that has been confirmed via a different approach [7]. This two 
dimensional analysis describes TE-polarized light in a cylindrical medium (where the electric- 
field vector E points orthogonal to the plane), and the electric field magnitude fulfills the 
Helmholtz equation. Leonhardt and Philbin have also demonstrated the analogous ideality of a 
novel impedance-matched spherical MFE for perfect focusing of electromagnetic waves in three 
dimensions [8]. 

In the two-dimensional case, the perfect focusing of the MFE in [6] assures that the medium 
will perfectly transport an outward (monopole) Helmholtz wave field, one generated by a point 
source, towards an "infinitely-well localized drain" [6] (one that we will call "perfect point 
drain") located at the desired image point. Note that the perfect point drain must be such that it 
totally absorbs all incident radiation, with no reflection or scattering by it. Note also that the 
field around the drain asymptotically coincides with an inward (monopole) wave. We will refer 
here to such a wave as "Leonhardt' s forward wave". 

Even though the physical significance of a point source as a limiting case seems to be well 
accepted, that of a passive perfect point drain has been considered very controversial [9] 
[10][1 1][12]. In reference [6], the drain was not physically described, but only considered as a 
mathematical entity, leaving no clues as to how to simulate that drain in software. Particularly, 
an analysis of such a drain, one located at a position different from the image point, would help 
to prove the super-resolution, which could not be done with the information in reference [6]. 



Recently, however, a candidate for perfect drain has been proposed for a microwave-frequency 
MFE [13], wherein a two-dimensional MFE medium has been assembled as a planar waveguide 
with concentric layers of copper circuit board forming the desired refractive index profile of the 
MFE. Also, both source and drain have been built as identical coaxial probes, one to introduce 
power into the planar waveguide and the other to extract it. The coaxial probe was intended to 
act as the perfect sink and is completely passive and loaded with the characteristic impedance of 
the coaxial cable [13]. There is no theoretical proof that such a coaxial probe will actually 
behave as a perfect drain, but a practical proof is claimed by comparing (Figure 4 in [13]) the 
measured electric field distribution and the analytical expression of Leonhardt's forward wave 
([6], reviewed here in Section 1.1). Since the measured and analytical values differ significantly 
(deviations attributed primarily to imperfections in the probes), we think a more detailed 
analysis is called for in order to clarify whether or not that receiving coaxial probe is acting as 
the perfect drain as defined here. Nevertheless, this clarification does not affect the main 
conclusion in [13]: such a probe certainly leads to super-resolution in the conditions of that 
measurement. 

In this paper we present a different realization of a passive perfect drain, obtained theoretically 
from the analytical equations. It consists in the introduction of a certain non-magnetic material 
inside a circle of radius R enclosing the image point (although not centered upon it, as discussed 
later). As will be proved in Section 2, the complex permittivity of that non-magnetic material 
will be such that the field outside that circle coincides exactly with Leonhardt's forward wave 
and the incoming power will be fully absorbed inside it. In general, this perfect drain can have a 
finite radius R; the perfect point drain is just the limit case R^>0. 



1.1 Leonhardt's forward wave 

The strength of the TE monochromatic field E(x,y)e'" a z in a region without sources or drains 
fulfils the 2D-Helmholtz equation: 

AE + n 2 k*E = (1) 
where k = a)^/u s a and n is the MFE refractive index distribution given by: 
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where p = ^jx 2 +y 2 . Leonhardt's forward wave is a particular family of solutions of Eq. (1), 
given by Eq. (12) in [6]: 

E= p J£ti^Hl (3) 

4 sin(v;r) 

where P v is the Legendre function of the first kind, 
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Here z=x+iy is the complex notation of the point (x,y), and x is an arbitrary real number. 
Without lack of generality we have located the point source described in [5] at the object point 
(x , 0). Note that -1<^<1 and by the divergence of PjiO when E is infinite at \Q=\. A 

wave according to Eq. (3) is generated by the point source located at (x ,0), and it propagates 
towards the perfect point drain at the image point (-l/x , 0). The time evolution of the field, 
Re(E(x,y)e~"°') is shown in the associated media files Mediala and Medialb for x =-2, £=15. 
This evolution clearly includes the one-directional propagation of the wave from source to 
drain, with no reflection or scattering by it. 

Note that the magnetic field H(x,y) of Leonhardt's forward wave can be easily computed from 
the field E, Eq.(3), and Eq. (5), as: 
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1 . 2 Current of the source 

Eq. (1) is only valid in a region without sources or drains, in our case the full plane except for 
the points (x , 0) and (-l/x , 0). The equation valid for the full plane must include the Dirac 
delta at those points. The amplitude of Eq. (3) was specifically selected in [6] to make it behave 
as a Green function, i.e., so the weights of those Dirac deltas have unit moduli. It can be written 
as [6][7][8]: 

AE + n 2 klE = -S(x — x , y) — c' m S(x — l/x ,y) (7) 

The right hand side of Eq. (7) can be identified (from Maxwell equations) as -z'ft>//oJ(x,y)z, 
where J(x,y) is the current density. Consequently the electric currents through source and drain 
for the specific field amplitude of Eq. (3) are: 
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1.3 Alternative expression for the forward wave 

There is an alternative way to express Eq. (3) (Eq. 12 of [6]), which uses the Legendre function 
of the second kind Q v . We consider here the branch of g„that is real-valued when lm(<^)=0 and 
|^1<1 (another complex -valued branch Q v is claimed to be considered in [6]). Taking into 
account, from [14] (Eq. (15) in p. 144), that: 

P v (-0 = cos(xv)P v (C)sm(xv)-Q v (<;) (9) 

n 

We already know that the Eq. (3)can be alternatively rewritten as: 
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Let us designate the factor inside the parenthesis of (10) as: 

F v (C) = P v (C) + i-Q v (C) (11) 
n 

Interesting asymptotic expressions for function F v are set forth in Appendix 1, which gives 
comprehensive discussion of Leonhardt's forward wave. Additionally, this alternative 
expression will greatly simplify some calculations in the perfect sink design of the next section. 



2 Perfect sink 



The theoretically ideal point drain of Leonhardt is located at the point (-1/xo, 0). We will design 
first a finite-area perfect drain, which will comprise a convex region surrounding that point 
(-1/xo, 0), filled with an inhomogeneous, isotropic, non-magnetic material with complex 
dielectric permittivity (thus, absorptive), such that the field outside that region coincides exactly 
with Leonhardt' s forward wave. 

The equation satisfied by the field in that region will then be the homogeneous Helmholtz 
equation. Since the incident wave fields E and H are known, the necessary continuity of E and 
H on the drain boundary will be forced by particularizing on the boundary the values of E and 
H in Eq. (3) and Eq.(6), respectively. 



2. 1 Selection of the boundary 

From Eq. (5), it is easy to confirm that the line £=Q =const is a circle with its centre at the point 
(x c , 0) and its radius R given by: 

2x R= -^c^ + 4)_ 

Xc l-x 2 -(l + x 2 )C d l-x 2 -(l + x % ° 2) 



We select the drain region as that containing the points fulfilling £>Q, with Q fulfilling: 

C^T-TJ (13) 
l + |x | 

Condition Eq. (13) is required in order to guarantee that the drain area is finite and encloses the 
image point (-l/x , 0), as can be deduced from Eq. (12). The ^constant circumferences are 
shown in red in Fig.l. The ^constant circumferences are also i?=constant lines (see Eq.(3)), 
and coincide with the GO wavefronts. This property is fulfilled by a more general class of 
inhomogeneous media, a class analyzed in [7]. The blue lines in Fig.l, which we can call 
<5=constant, are also circumferences, and coincide with the GO rays, and are the Poynting vector 
lines of Leonhardt's forward wave. Both families of lines coincide with coordinate grid lines of 
the bipolar orthogonal coordinate system (see, for instance, [7]). 

From Eq. (12) we see that selecting Q close enough to 1 will make the radius R as small as 
desired. At the limit we obtain the point drain x c — >■ -l/x ). 




Fig. 1 . The source is a line of current perpendicular to the x-y plane and placed at the 
point (x , 0), which defines the two-dimensional point source. The boundary of the 
finite-area passive drain is a cylinder with base Q> The red lines (^constant) and 
the blue lines (^constant) define a cylindrical bipolar coordinate system. 



2.2 Inhomogeneous complex refractive index of the drain 

Inside the drain (£> Q), we select // =// and the refractive index has the following form: 

"' = TV <14) 

where Ed is a complex constant and with Im^^O to be calculated later (in Section 3). Then, the 
homogeneous Helmholtz equation in the drain is: 

AE + n 2 XE = (15) 

Designate k d = a>^ju s s d (which is complex). Using the expression for the refractive index n of 
the MFE (Eq. (2)), we find that the selection made with Eq. (14) fulfills that n d k =nk d , so Eq. 
(15) can also be written as: 

AE + n 2 kjE = (16) 

This equation is identical to the Helmholtz equation of the MFE, Eq.(l), after substituting the 
real wave number k by k d (still to be calculated). 

2. 3 Ordinary differential equation of the drain 

One of the boundary conditions on the line Q is the continuity of field E=Ez. As said before, 
E is constant on that boundary surface, so it would be particularly interesting to express Eq. (16) 
in the bipolar coordinate system £-8. That was already done in Section 3.1 of reference [7]. As 
shown there, the expression of Eq. (16) for solutions depending only on <^is the same as that Eq. 
(9) in [6]. Using the change of variables <^=(r 2 -l)/(r 2 +l), this equation is: 
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Its general solution [14] can be written as: 

E = AP V ,(C)+BQ,{C) 

where V is given by 



(18) 



2 



(19) 



The three constants A, B and v are fixed by three conditions. Two of them are given by the 
continuity of the E and H fields at the boundary, and the third is that the field E must be 
bounded (i.e., it cannot diverge), since the Helmholtz Eq. (15) is homogeneous. 

Consider first the third condition. From the properties of the Legendre functions [14], we know 
that the function Q v {£) diverges when if lm(v')^0 and P v {£) does not (P v {\)=\). 

Therefore, this boundary condition imposes B=0, which means that the field inside the drain 
region (<^> <^ d ) has the form: 

E = AP v ,(£) (20) 



A and v are calculated forcing the other two conditions i.e., the continuity of E and H at the 
boundary. E and H outside the drain are taken from the solution in the absence of reversed wave 
(so the drain is reflection-less) i.e., by using Eq. (10) and Eq. (6)). 
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Dividing both equations, we obtain: 
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(22) 



which is an expression where only v is unknown. Once this is solved we can calculate e d with 
Eq. (19) : 



e d = 



J 



v'(v'+l) 



(23) 



3. Examples 



Fig. 2 shows s d as function of frequency for two different drain radii and Fig. 3 shows it as a 
function of the drain radius for different frequencies. In all the examples we shall assume that p 



in Eq (2) and (14) is given in cm and frequency in GHz. Media2a.mov and Media2b.mov show 
the time evolution of the field, Re(E(x,y)e'" a> ) for a drain radius i?=0.2cm, Xo=-2, k=\5, This 
should be compared with the case of the perfect point drain shown in the Introduction (Media2a 
and Media2b). Media3 shows the time evolution of the field inside the drain for radius 0.2 cm. 
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Fig. 2. Dependence of £a (Eq. (14) with frequency. Dark and bright blue curves 
correspond to the real and imaginary part for a drain radius 0.2 cm centered at 
x c =0.5159 cm, y c ~0. Red and orange curves correspond to the real and imaginary parts 
for a drain radius 0.001 cm centered at x c =0.5 cm, j c =0. The point source is located at 
x =-2 cm, yo=0. Frequency in GHz. 
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Fig. 3. Dependence of Ej (Eq. (14) with drain radius. Dark and bright blue curves 
correspond to the real and imaginary part for a frequency 200GHz. Red and orange 
curves to real and imaginary part for a frequency 10GHz. 



4. Looking inside the drain. 

4.1. Electric field and current inside the drain. 



The conductivity a of the media inside the drain can be calculated as a function of £j from its 
definition: 



CO 

a - — Im 
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(24) 



Using Eq. (23) we can calculate eras a function of v. The electric field inside the drain is given 
by Eq. (20). Consequently the current density J=<tE in the drain is: 



J(x,y) = - ^Im[£jW)z (25) 

0-+P ) 



Fig. 4 shows the modulus of the current density inside the drain for a drain radius=0.2 cm and a 
total current in the source of 1mA when the drain center is at the point (0.519, 0) cm. 




Fig. 4. Representation of the modulus of the current density (mA/cm 2 ) inside the drain. 
Drain radius=0.2 cm, drain center at (0.519,0) cm. Total current in the source =lmA. 

4.2 Power absorption. 

The power emitted by the source P can be obtained by integrating the Poynting vector over a 
surface enclosing the source. The MFE is a lossless system, so this power has to equal both the 
total power entering the drain and the power it absorbs. This integration has been made as 
shown in Eq. (26). This surface has cylindrical symmetry along the z-axis, so for the sake of 
simplicity we take a surface whose projection on the x-y plane is a ^=constant curve, let's say 
C = Cs- Since the Poynting vector has no z-dependence the surface integral is reduced to a line 
integral along Q=£ s . 



P = -Re[jExH* 
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Alternatively, the same power can be calculated as the power absorbed in the volume of the 
drain, which is: 



Hi 



J-EdV 



J cr\AP v ,(0\ 2 dS 



= -Re 
2 



dPAO 




* 


















dx y 



dl 



(27) 



Fig. 5 shows the power P as function of the frequency for 1mA total current in the point source. 
Note this implies an amplitude value for Leonhardt's forward wave, correspondingly scaled 



from Eq. (3) and (8). The power P depends on the source total current and frequency, but 
obviously it does not depend on the drain radius. 
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Fig. 5 Power emitted by the source and absorbed by the drain as function of frequency. 
Total current in the source is 1 mA. 

5. Conclusions. 

We have found that the perfect drain concept can be modeled as a dissipative region the 
diameter of which tends to zero. The Leonhardt's forward wave solution inside the MFE can not 
be obtained for any material used as a drain. We have calculated the complex permittivity of the 
non-magnetic material forming the perfect drain. When the size of the drain tends to zero, the 
drain tends to the perfect drain concept Introduced by Leonhardt. Without such a perfect drain, 
there will be two waves Fj[Q and RXO inside the MFE (forward and backward respectively, 
and described in Eq.(3 1)) as response to the point source. 

This concept of the drain as a small dissipative region can be easily included in electromagnetic 
modeling software. Both characteristics are desirable to analyze and experimentally verify the 
super-resolution properties of the MFE, i.e. to check that the power P changes drastically when 
the source (or the drain) moves to a neighbor point located at a distance much smaller than the 
wavelength. 

Appendix. Asymptotic expression when v »1, and the Backward wave. 

Using Eq. (1) and Eq. (2) in page 162 of [14], functions PJ^Q and QJ^Q, £ = cosd, for s 
<9<?r-e(e>0) can be approximated as: 

P v (cos 0) = F(V) , /— (cos (( v + 1 / 2) 0) + O ( v- 1 )} 
vV ' r(v + 3/2)Vsin^ vv i ) \ )\ 

, (28) 

Q v (cos 0) = F(V) J— (- sin (( v + 1 / 2) 0) + O (v~ x )) 
* vK ' I> + 3/2)V2sin^ U ) ) \ )S 



Thus, using Eq.(l 1) we have: 



F (cos 0) = F(v) J— L<»v*» + o (v- 1 )} (29) 
vK } I> + 3/2)V;z-sin^ 1 '/ V ' 

which is clearly identified as a wave propagating towards decreasing values of 9 (remember the 



factor e'"*), that is, increasing £ from source to drain. This expression is accurate for 
s<9<n-e (s >0) (that is, not too close to the point source and point drain). Asymptotes to Eq. 
(28) and (29) suggest to define the function: 

R v (C) = P v (C)--jQ v (C) (30) 

71 

to describe a field propagating from image point to object point. In analogy to canonical 
solutions of the Helmholtz Equation (as plan waves, cylindrical waves or spherical waves in 
free-space), it seems more natural for this MFE problem to define the general solution of Eq. 
(17) as a superposition of functions F v and R v , that is, to rewrite the field E as : 

E = CF v (C) + DR v (C) (31) 

When there is a perfect drain, outside it only the Leonhardt's forward wave F v exists, which 
implies that D in Eq. (31) is null. Inside the drain, the condition that the field E must be 
bounded leads to Eq. (20), which in terms of F v - and R v implies the condition C=D. This means 
that inside the perfect drain there is a standing wave, which is necessary to avoid the singularity 
at the image point (<^=1). This is analogous to the superposition of converging and diverging 
cylindrical waves, of equal amplitude, in free-space as described by the Hankel functions 

H ( l) (k p) and H q 2) (k p) , which results in the bounded Bessel function Jo(k p). 

The forward and backward wave defined here for the MFE lens are intimately related to the 
retarded and advanced field defined in [12] and [13] for the MFE mirror. Using the formulation 
in [12] and [13], a wave bounded at the image point for the MFE lens can be written: 

f? — „ fe ( v (*b)-v(-*<>)) 77 

j wr(v(*o)-v(-*o)) 1 ' 

where the positive and negative wave-numbers are defined as 

-l±Jl + 4£ 2 

K±*b) = \ (33) 

Besides the complexity of expression Eq. (32), it can be easily seen that it is (up to a 
multiplicative constant) simply equal to P v (which is the bounded expression used above in Eq. 
(20)). This is obtained by direct computation, after considering that v(-k )= -v(k )-l, and then 
using Eq. (7) and (16) in page 144 of [14]: 
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